Ultra-High Energy Cosmic Rays in a Structured and Magnetized Universe 



m 
o 
o 

(N 



Giinter SigP, Francesco Miniati'', Torsten A. EnBlin'' 
GReCO, Institut d'Astrophysique de Paris, C.N.R.S., 98 bis boulevard Arago, F-75014 Paris, France 
* Max-Planck Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, 8574-1 Garching, Germany 

We simulate propagation of cosmic ray nucleons above 10^^ eV in scenarios where both the source 
distribution and magnetic fields within about 50 Mpc from us are obtained from an unconstrained 
large scale structure simulation. We find that consistency of predicted sky distributions with current 
data above 4 x 10^^ eV requires magnetic fields of ~ OAfiG in our immediate environment, and a 
nearby source density of ~ 10~^ — 10"'^ Mpc""^. Radio galaxies could provide the required sources, 
but only if both high and low-luminosity radio galaxies are very efficient cosmic ray accelerators. 
Moreover, at ~ 10^® eV an additional isotropic fiux component, presumably of cosmological origin, 
should dominate over the local flux component by about a factor three in order to explain the 
observed isotropy. This argues against the scenario in which local astrophysical sources of cosmic rays 
above ~ 10^^ eV reside in strongly magnetized {B ~ O.l/xG) and structured intergalactic medium. 
Finally we discuss how future large scale full-sky detectors such as the Pierre Auger project will 
allow to put much more stringent constraints on source and magnetic field distributions. 

PACS numbers: 98.70.Sa, 13.85.Tp, 98.65.Dx, 98.54.Cm 
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I. INTRODUCTION 



Over the last few years the detection of several giant 
air showers, either throughground based detectors 0, Q 
or fluorescence telescopes 0,131: has confirmed the arrival 
of ultra high energy cosmic-rays (UHECRs) with energies 
up to a few hundred EeV (1 EeV = 10^* eV). Their exis- 
tence poses a serious challenge and is currently subject of 
much theoretical research as well as experimental efforts 
(for recent reviews see 0,0,0)- 

The problems encountered in trying to explain UHE- 
CRs in terms of "bottom-up" acceleration mechanisms 
have been well-documented in a number of studies (e.g., 
Refs. 0, 01). In summary, apart from the specific 
energy draining interactions in the source the maximal 
UHECR energy is limited by the product of the accelera- 
tor size and the strength of the magnetic field. According 
to this criterion it turns out that it is very hard to acceler- 
ate protons and heavy nuclei up to the observed energies, 
even for the most powerful astrophysical objects such as 
radio galaxies and active galactic nuclei. 

In addition, nucleons above ~ 70 EeV suffer heavy en- 
ergy losses due to photo-pion production on the cosmic 
microwave background (CMB) — the Greisen-Zatsepin- 
Kuzmin (GZK) effect |il| — which limits the distance 
to possible sources to less than ~ 100 Mpc 0. Heavy 
nuclei at these energies a re p hoto-disintegrated in the 
CMB within a few Mpc [l^l- Unless the sources are 
strongly clustered in our local cosmic environment, a 
drop, often called the "GZK cut-off" in the spectrum 
above ~ 70 EeV is therefore expected , even if the in- 
jection spectra extend to much higher energies. However, 
the existence of the latter is not established yet from the 
observations In fact, whereas a cut-off seems con- 

sistent with the few events above 10^° eV recorded by 
the fluorescence detector HiRes 0, it is not compatible 
with the 8 events (also above 10^° eV) measured by the 
AGASA ground array 0|. The solution of this problem 



may have to await the completion of the Pierre Auger 
project jl6i| which will combine the two complementary 
detection techniques adopted by the aforementioned ex- 
periments. 

Adding to the problem, there are no obvious astronom- 
ical counterparts to the detected UHECR events within 
~ 100 Mpc of the Earth 0, At the same time, 

no significant large-scale anisotropy has been observed 
in UHECR arrival directions above ~ 10^^ eV, whereas 
there are strong hints for small-scale clustering: The 
AGASA experiment has observed five doublets and one 
triplet within 2.5° out of a total of 57 events detected 
above 40 EeV 01 ■ When combined with three other 
ground array experiments, these numbers increase to at 
least eight doublets and two triplets within 4° T^. This 
clustering has a chance probability of less than 1% in the 
case of an isotropic distribution. 

Independent of the specific UHECR production mecha- 
nism, there are currently two possible explanations of the 
experimental findings described above: The first assumes 
very weak intergalactic magnetic fields capable of deflect- 
ing UHECRs only up to a few degrees, or neutral pri- 
maries. In this case the apparent isotropy would indicate 
that many sources contribute to the observed flux and 
most of these sources would be at cosmological distances 
because the local source distribution is in general too 
anisotropic to be consistent with the observed UHECR 
isotropy. This would also explain the absence of nearby 
counterparts and a subset of especially powerful sources 
would explain the small-scale clustering Indeed, it 
has been argued that UHECR arrival directions corre- 
late with the positions of BL Lacertae objects, suggesting 
these as sources accelerating protons j23| , although there 
seems to be disagreement about this in the literature . 
Furthermore, some of these objects may be too far away 
to be consistent with the GZK effect, which would req uire 
new physics such as Lorentz symmetry violations |23 | . In 
contrast, correlations with compact radio quasars have 
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not been found |23|- If correlations with astrophysical 
objects are confirmed, this would strongly suggest small 
deflection or neutral primary particles. Whatever the 
sources are in this scenario, for small deflection one can 
in principle constrain the characteristics of the magnetic 
fields along the line of sight and the source properties 
by analysing arrival times, directions, and energies of ob- 
served small-scale multi-plets 01 • Also, in the small de- 
flection scenario the experimental conflrmation of a GZK 
cutoff is expected. 

However, the assumption of weak intergalactic mag - 
netic fields seems at odds with several observations [25j . 
Most remarkable are the detections of Faraday rotation 
measures which seem to indicate field strengths at the 
/iG level within the inner region central Mpc) of 
galaxy clusters l^^. In addition, the recent mounting 
evidence for diffuse radio-synchrotron emission in numer- 
ous galaxy clusters W?\ and in a few cases of filaments 
[28I I29I I , seems to suggest the presence of magnetic fields 
as strong as 0.1-1.0/iG at the relatively low density out- 
skirts of collapsed cosmological structures. In fact, extra- 
galactic magnetic fields (EGMF) as strong as ~ 1/iG in 
sheets and filaments of the large scale galaxy distribution, 
such as in our Local Supercluster, are compatible with 
existing upper limits on Faraday rotation [23, Isol l3ll | . It 
is also possible that fossil cocoons of former radio galax- 
ies, so called radio ghosts, contribute to the isotropiza- 
tion of UHECR arrival directions [32|. Thus, relatively 
strong magnetic fields seem to be ubiquitous in inter- 
galactic space, although their theoretical understanding 
is still hmited |3|. 

Such observational evidence motivates a second, more 
realistic scenario, which takes into account the exis- 
tence of strong (S ^ 0.1 — 1 /iG) intergalactic magnetic 
fields correlated with the large scale structure. In this 
case magnetic deflection of charged primaries would be 
considerable even at the highest energies and the ob- 
served UHEGR flux could be dominated by relatively 
few sources within about 100 Mpc. Here, large scale 
isotropy could be explained by considerable angular de- 
flection leading to diffusion up to almost the highest en- 
ergies and the small scale clustering could be due to mag- 
netic lensing [s^]. The locations of clusters of events of 
different energies would in this case coincide with the 
crossing points of the caustics for these energies where 
fluxes are enhanced. 

In the present paper we take this second point of view 
and investigate in some detail the effects of propagation 
of UHECRs, assumed to be dominantly nuclcons, in a 
magnetized large scale structure matter distribution com- 
puted according to a numerical cosmological simulation. 

Early investigations of this scenario have been carried 
out in Refs. 35, 36, 37, 38, 39j, assuming that sources and 
magnetic fields follow a pancake profile of scale height 
~ 3 Mpc and scale length ~ 20 Mpc, the magnetic field 
having a power law spectrum at length scales below 
~ 1 Mpc. UHECR propagation was computed through 
a numerical code that accounts for magnetically induced 



defiections and all relevant energy losses [35*, "3^, "sT] . The 
cases of a single source js^ |^ , as well as continuous ^3^ 
and discrete source distributions 39] have been investi- 
gated. The above studies led to the result that the multi- 
pole moments and autocorrelation functions of the arrival 
directions best fit the AG AS A data for a number ^10 
sources in the Local Supercluster, assumed to emit con- 
tinuously, and a maximal field strength of ~ 0.3/iG 3^. 

Ideally, however, it would be desirable to study the 
propagation of UHECRs based on distributions of both 
potential sources and observed magnetic field properties. 
However, up to now, only catalogs of candidate sources 
have been available. Magnetic fields, on the other hand, 
have been approximated in a number of fashions: as neg- 
ligible ^3 ) as uniform [ij , or as organized in spatial cells 
with a given coherence length and a strength depending 
as a power law on the local density p^. 

In the present paper we attempt to go beyond some 
of the above limitations by computing for the first time 
the propagation of the UHECRs in a magnetized cosmo- 
logical environment computed through numerical simu- 
lations. We carry out a fully cosmological simulation 
of large scale structure formation which, in addition to 
dark matter and baryonic gas, follows the evolution of 
a passive magnetic field. This approach is motivated by 
the fact that /iG magnetic fields are mostly negligible 
for the purpose of the dynamics of the large scale cos- 
mic fiows (hence their passive character). In addition, 
and basically for the same reason, the structure of mag- 
netic fields on scales of interest for UHEGR propagation 
(~ 100 kpc) is mostly determined by the hydrodynamic 
flow. This is conflrmed by the fact that in these sim- 
ulations, the magnetic field looses memory of its initial 
conditions, soon after the formation of structures begins. 
Finally, the statistical properties of cosmological struc- 
ture in the universe are rather homogeneous. Therefore, 
the simulated matter structure and magnetic field distri- 
butions should provide a realistic scenario for studying 
the statistical properties of UHEGR source distributions 
and propagation in a cosmic environment. In the present 
study we assume the sources to follow the baryon den- 
sity. Furthermore, the observer is supposed to be in re- 
gions of the simulated matter distribution which contain 
structures of the same size and baryonic gas tempera- 
ture as our local neighborhood. This should provide a 
suitable environment to simulate the arrival of UHECRs 
from extragalactic distances and the effects of local mag- 
netic fields of various strengths. 

In the future such studies can be further improved by 
computing constrained simulations that reproduce in de- 
tail the observed matter distribution of the local uni- 
verse. Such a simulation has been used for the case of 
radio ghosts in Ref. _32J where, however, the magnetic 
fields were not followed but were rather assumed to scale 
with the gas density. Constrained simulations including 
magnetic fields are relevant for predicting quantitative 
features such as location of clustered events, phases of 
anisotropics etc. and will be used in a following study. 
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We point out, however, that for the reasons given above, 
effects of the magnetic field and source distributions in 
the local universe should essentially be captured by the 
present approach at least up to "cosmic variance" . The 
latter represents variations due to different source and 
observer locations and will be estimated in our simula- 
tions. 

We also restrict ourselves to UHECR nucleons, and 
we neglect the Galactic contribution to the deflection of 
UHECR nucleons since typical proton deflection angles 
in galactic magnetic fields of several are ^ 10° above 
4 X 10^^ eV ^31 J ^-iid thus in general are small compared 
to extra-galactic deflection in the scenarios studied in the 
present paper. 

The simulation is described in more detail in the next 
section. There we also describe the general features of 
our method and define the statistical quantities used for 
comparison with the data. In Sect. 3 we present results 
and we conclude in Sect. 4. 



II. MOTIVATION AND OUTLINE OF THE 
NUMERICAL MODEL 

A. Magnetic Deflection 

Contrary to the case of electrons, for charged hadrons 
deflection is more important than synchrotron loss in 
the EGMF. To get an impression of typical deflection 
angles one can characterize the EGMF by its r.m.s. 
strength B and a coherence length Ic- If we neglect en- 
ergy loss processes for the moment, then the r.m.s. de- 
flection angle over a distance r ^ in such a field is 
9{E,r) ~ (2r/c/9)^/^/rL 0' where the Larmor radius 
of a particle of charge Ze and energy E isr^ — E/ [ZeB). 
In numbers this reads 

for r ^ Ic- This expression makes it immediately obvi- 
ous why a magnetized Local Supercluster with fields of 
fractions of micro Gauss prevents a direct assignment of 
sources in the arrival directions of observed UHECRs; 
the deflection expected is many tens of degrees even at 
the highest energies. This goes along with a time delay 

T{E,r) ~ re{E,df/A (2) 

Vio2"evy vioMpcy 

'''' 

which may be millions of years. A source visible in UHE- 
CRs today could therefore be optically invisible since 



many models involving, for example, active galaxies as 
UHECR accelerators, predict variability on much shorter 
time scales. 



B. Numerical Simulation of the Large Scale 
Structure 

The formation and evolution of the large scale struc- 
ture is computed by means of an Eulerian, grid based 
Total- Variation-Diminishing hydro-|-N-body code [45l |. 
We adopt a canonical, flat ACDM cosmological model 
with a total mass density = 0.3 and a vacuum en- 
ergy density JIa = 1 — — 0.7. We assume a nor- 
malized Hubble constant h^i = Ho/67 km s^^ Mpc^^ 
~ 1 and a baryonic mass density, Qb — 0.04. The sim- 
ulation is started at redshift z ~ 60 with initial density 
perturbations generated as a Gaussian random field and 
characterized by a power spectrum with a spectral index 
Us = I and "cluster-normalization" erg = 0.9. 

We adopt a computational box size of 50 Hq^^ Mpc. In 
this box the dark matter component is described by 256'^ 
particles whereas the gas component is evolved on a co- 
moving grid of 512^ zones. Thus each numerical cell 
measures about 100 hg^^ kpc (comoving) and each dark 
matter particle corresponds to 2 x 10^ Hq^ Mq. Besides 
the box and dark matter particle sizes the cosmological 
simulation is the same as that presented in Ref . • 

The magnetic field is followed as a passive quantity, 
that is magnetic forces are neglected. This is consistent 
with the strength of observed magnetic fields in most dif- 
fuse extragalactic environments. Basically we solve the 
induction equation with the velocity field provided by the 
simulated flow |^ and the initial magnetic field seeds 
generated by the Biermann battery mechanism. How- 
ever, as already pointed out, the initial conditions are not 
important as the topological properties of the magnetic 
field are determined by the subsequent evolution of the 
large scale fiow. This is responsible for its amplification 
through gas compression and shear flows. Thus, at the 
end of the simulation, the relative strength of the mag- 
netic field in different regions is determined by the hydro- 
dynamic properties of the flow. While the simulation out- 
come regarding the relative magnetic field strength and 
topology distribution are obviously retained, the overall 
normalization is chosen in order to reproduce the fields 
of several micro Gauss observed in the regions of largest 
density, namely galaxy cluster cores. Fig. Q illustrates 
an example of the simulated magnetic pressure (top) 
and baryonic density (bottom) distributions. The figure 
shows two-dimensional cuts corresponding to a depth of 
100/igy^kpc. The color images are in log scale and, for 
visualization purposes, span a dynamic range of 3 and 
6 orders of magnitude for magnetic pressure and bary- 
onic density respectively. The magnetic field is partic- 
ularly strong in both postshock regions and inside rel- 
atively large structures where it has been compressed 
and stretched. Apparently, its distribution is less con- 
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FIG. 1: Log-scale two-dimensional cut through magnetic 
pressure (top) and baryon density (bottom). The image is 
50 hg^ Mpc on each side and 100 h^^ kpc deep. The small 
white dot in the bottom panel indicates the location of the ob- 
server. For visualization purposes we adopt a dynamic range 
of 3 and 6 orders of magnitude for the magnetic pressure and 
baryon density, respectively. 



centrated than the baryonic density, resembhng in this 
respect that of the thermal pressure (not shown). 



C. Simulated UHECR Experiments 

To simulate the propagation and arrival of UHECRs in 
the computational box we need to choose: (a) the loca- 
tion of the observer and (b) the source distribution. As 



anticipated in the introduction, the location of the ob- 
server is identified as a region whose general features in 
terms of scale, mass and temperature, resemble those of 
the local universe. That means a small group of galaxies 
characterized by a gas temperature of order of a frac- 
tion of a keV. There are several such structures in a 
50 hgj Mpc box such as the one employed here. In the 
neighborhood of the one we selected as the observer lo- 
cation, we also find a larger group of galaxies with tem- 
perature of a few keV. In order to orient the simulation 
box with respect to the observed sky, the latter object, 
located at a distance of ~ 34 Mpc, is arbitrarily asso- 
ciated with the Virgo cluster. This reference frame al- 
lows us to define a celestial system of coordinates {a, S) 
which describes the arrival direction of events recorded 
by our virtual observer. It will be useful in the next sec- 
tion where the arrival direction probability distribution 
is constructed. The above setting is sufficient for the 
current purpose of investigating the effects on the prop- 
agation of UHECRs of realistic, topologically structured 
magnetic fields of various strengths. 

We then chose randomly a certain total number Ng of 
sources in the box, corresponding to an average source 
density 8 x 10~^/i57 TV^ Mpc~'^, with probability propor- 
tional to the local baryon density. In order to avoid intro- 
ducing too many free parameters, we further assume that 
all sources roughly emit the same power law spectrum of 
CRs extending up to ~ 10^^ eV, with roughly equal to- 
tal power. We also assume that neither total power nor 
the power law spectral index change significantly on the 
time scale of UHECR propagation. This can be up to a 
few Giga years for the magnetic fields considered here. 
Injected power and spectral index are then treated as 
parameters which can be fit to reproduce the observed 
spectrum, as will be seen below. 

For each such configuration many nucleon trajectories 
originating from the sources were computed numerically 
by solving the equation of motion for the Lorentz force 
and checking for pion production every fraction of a Mpc 
according to the total interaction rate with the CMB and, 
in case of an interaction, by randomly selecting the sec- 
ondary energies according to the differential cross section. 
Pair production by protons is treated as a continuous en- 
ergy loss process. 

A detection event was registered and its arrival direc- 
tions and energies recorded each time the trajectory of 
the propagating particle crossed a sphere of radius 1 Mpc 
around the observer. For each configuration this was 
done until 5000 events where reg istered. For more de- 
tails on this method see Refs. |35[l36ll37| . 



D. Data Processing 

For each realization of sources and observer, these 
events were used to construct arrival direction probability 
distributions, taking into account the solid- angle depen- 
dent exposure function for the respective experiment and 
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folding over the angular resolution. 

For the exposure function uj{d) we use the parameter- 
ization of Ref. ^3 which depends only on declination 
S, 

u>{5) oc cos(ao) cos(5) sin(am) + ctm sin(ao) sin(5) , 

[ if e > 1 

where = < tt if ^ < — 1 , (3) 

[ cos~^(^) otherwise 

. , , ^ _ cos(6'„) - sin(ao) sm{S) 

witli 4 = — — . 

cos(ao) cos(o) 

For the AGASA experiment aq = -35°, 0„i = 60°, and 
the angular resolution 2.4° are used. For a full-sky Pierre 
Auger type experiment we add the exposures for the 
Southern Auger site with gq — —35° and a putative sim- 
ilar Northern site with gq — 39°, and 6m = 60° in both 
cases, with an assumed angular resolution of ~ 1°. 

From the distributions obtained in this way typically 
1000 mock data sets consisting of N observed events were 
selected randomly. For each such mock data set or for the 
real data set we then obtained estimators for the spher- 
ical harmonic coefficients C{1) and the autocorrelation 
function N{9). The estimator for C{1) is defined as 

1 1 ' / " 1 \' 

m= — l \i=l / 

where uji is the total experimental exposure at arrival 
direction u^, — X^Hi ^/^i is the sum of the weights 
and Yimiu'') is the real- valued spherical harmonics 
function taken at direction u*. The estimator for N{d) is 
defined as 

_ C J 1 if 9ij is in same bin as 1 
^ ' ~ ^10 otherwise j ' 

and S{6) is the solid angle size of the corresponding bin. 
In Eq. © the normalization factor C = ile/{N{N - 1)), 
with rig denoting the solid angle of the sky region where 
the experiment has non-vanishing exposure, is chosen 
such that an isotropic distribution corresponds to = 
1. 

The different mock data sets in the various realizations 
yield the statistical distributions of CiX) and N{6). One 
defines the average over all mock data sets and realiza- 
tions as well as two errors. The smaller error (shown to 
the left of the average in the figures below) is the statis- 
tical error, i.e. the fluctuations due to the finite number 
N of observed events, averaged over all realizations. The 
larger error (shown to the right of the average in the fig- 
ures below) is the "total error" , i.e. the statistical error 
plus the cosmic variance. Thus, the latter includes the 
fluctuations due to flnite number of events and the varia- 
tion between different realizations of observer and source 
positions. 

Given a set of observed and simulated events, after 
extracting some useful statistical quantities Si, namely 



Ci and N{9) deflned above, we define 

i AS- J ' 

where Si^data refers to Si obtained from the real data, 
and S'i.simu and A^i^simu are the average and standard 
deviations of the simulated data sets. This measure of 
deviation from the average prediction can be used to ob- 
tain an overall likelihood for the consistency of a given 
theoretical model with an observed data set by counting 
the fraction of simulated data sets with Xn larger than 
the one for the real data. 



III. RESULTS 

In the following we compare the results obtained 
for the simulated UHECR propagation experiments de- 
scribed above with the observational results. In accord 
with what was outlined in the previous section, the com- 
parison is based on the statistical properties of the sim- 
ulated and observed events, expressed in terms of the 
angular power spectrum and the autocorrelation func- 
tion of the UHECR arrival distributions. A summary 
of the simulations run is contained in Tab. ^ There, 
for comparison, simulations 2 and 5 were performed for 
an observer situated in a small void with weak ambient 
magnetic fields. 

TABLE I: List of UHECR propagation simulations. The 
columns contain the simulation number, the number of 
sources in the simulation box of (50Mpc h^^l)^, the magnetic 
field strength at the observer location, the best fit power law 
index in the injection spectrum E~" , and the overall like- 
lihoods of fits to the AGASA data above 4 x 10^^ eV for 
the multi-poles Eq. I^J with I < 10 and the auto-correlation 
Eq. for 6 < 10°, respectively. The likelihoods are com- 
puted for n = 4 in Eq. @ which leads to reasonable discrim- 
inative power. 
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1 


100 


1.3 X lO"'^ 


2.4 


0.13 


0.63 


2 


100 


8.2 X 10"" 


2.7 


0.098 


0.15 


3 


10 


1.3 X lO"'^ 


2.4 


0.12 


0.69 


4 


10 


2.7 X lO"'^ 


2.4 


0.071 


0.15 


5 


10 


8.2 X 10"" 


2.7 


0.011 


0.037 


6 


1 


1.3 X lO"'^ 


2.8 


0.074 


0.62 



We find that as long as the observer is surrounded by 
magnetic fields of about O.lfiG, Ng ^ 10 nearby sources, 
i.e. sources within the simulation box, are necessary 
to reproduce multi-poles and autocorrelations marginally 
consistent with present data, limited, we emphasize, to 
the Northern hemisphere only. However, consistency 
of large scale multi-poles is somewhat worse than for 
the spatially more extended EGMF assumed in previ- 
ous work |39j. In Figs. [21 and |31 we show as an example 
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N„i„( 40. < E/EeV < 9000.)= 57, £= 0.185 




FIG. 2: The angular power spectrum C(l) as a function of 
multi-pole I, obtained for the AGASA exposure function, see 
text, for N — 57 events observed above 40 EeV, sampled from 
12 simulated configurations of scenario 1 in Tab. n The dia- 
monds indicate the realization averages, and the left and right 
error bars represent the statistical and total (including cosmic 
variance due to different realizations) error, respectively, see 
text for explanations. The histogram represents the AGASA 
data. The overall likelihood significance is ~ 0.13 for n = 4 
and / < 10 in Eq. ®. 



40- < E/EeV < 9000,)= 57, £= 0.746 




6 8 10 12 
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FIG. 3: As Fig. |5| but for the angular correlation function 
N{6) as a function of angular distance 6, using a bin size 
of A6 — 1°. Note that an isotropic distribution would cor- 
respond to N{6) — 1. The overall likelihood significance is 
~ 0.63 for n = 4 and 6 < 10° in Eq. It is not significantly 
diff'erent for somewhat larger bin sizes A6 ~ 2°. 



the results for the case of Ng = 100 nearby sources, sce- 
nario 1 in Tab. corresponding to a source density of 
8 X lO^^hl^ Mpc"^ The overall likelihood for ?i = 4 
in Eq. © is ~ 0.13 and ~ 0.63 for the multi-poles and 
autocorrelations shown, respectively. Also Fig. ^ shows 
that, for UHECR sources characterized by a proton in- 
jection spectrum roughly as oc E~^ '^ and extending up 



FIG. 4; Predicted spectrum observable by AGASA for 
the scenario 1 in Tab. |l] for which multipoles and auto- 
correlations were shown in Figs. |5| and |3 averaged over 12 
realizations, as compared to the AGASA (dots) and HiRes-I 
(stars) data. 



to ~ 10^^ eV, the observed spectrum at sub-GZK ener- 
gies is well reproduced. In addition, above GZK energies 
the spectral slope is predicted to be somewhere between 
the AGASA and HiRes observations, see Fig. 2| Nor- 
malizing to the observed flux results in a UHECR power 
of 5 X lO^^ergs^^ per source to be continuously emitted 
above 10^^ eV. 

The situation for Ns = 10 nearby sources does not lead 
to significantly different likelihoods, see scenarios 3 and 4 
in Tab. ^ However, the case of just one source is clearly 
disfavored in terms of the multi-poles, see scenario 6 in 
Tab.Q] This confirms similar findings in earlier work '36l. 

If the observer is in a region of EGMF strength much 
smaller than ~ 0.1/iG, as in scenario 2 of Tab. ^ for 
Ns ^ 100 nearby sources the predicted UHECR sky dis- 
tribution reflects the highly structured large scale galaxy 
distribution, smeared out only by the fields surround- 
ing the sources. This becomes obvious from Fig. [51 
which shows that UHECR arrival directions are much 
less isotropic in this case than if the observer is immersed 
in fields B ~ O.l/LtG. 

Nevertheless, the overall likelihood significance for 
multi-poles up to / = 10 is ~ 0.1, and thus not sig- 
nificantly worse than for the strong observer field case 
of Fig. [3 Therefore, the number of events observed 
by AGASA above 40 EeV is insufficient to distinguish 
this low observer field case from the strong observer field 
case based on anisotropy alone. However, as can be seen 
from Fig. |H| the low observer field case results in auto- 
correlations at angles 0^3° much larger than observed 
by AGASA. This is because strong magnetic fields at the 
observer position cause enough UHECR diffusion that 
their large-scale auto-correlations are significantly sup- 
pressed, as in Fig. |21 However, for fields considerably 
larger than 0.1/iG the auto-correlations tend to become 
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FIG. 5: Illustration of the influence of magnetic fields sur- 
rounding the observer on UHECR arrival direction distri- 
butions above 40 EeV in terrestrial coordinates. The up- 
per panel is for scenario 1 (observer surrounded by relatively 
strong magnetic flelds) , and the lower panel for scenario 2 (ob- 
server surrounded by negligible magnetic fields) from Tab.|T] 
averaged over all 12 and 10 realizations of 5000 trajectories 
each, respectively, thus corresponding to an effective number 
of sources of ~ 1000. The color scale represents the integral 
flux per solid angle. The pixel size is 1° and the image has 
been convolved to an angular resolution of 2.4° correspond- 
ing to the approximate AGASA angular resolution. The flUed 
sphere represents the position of the Virgo-like cluster. 



too strong again, see scenario 4 in Tab. U probably due 
to increased magnetic lensing. 

We also find that sources outside our Local Superclus- 
ter do not contribute significantly to the observable flux if 
the observer is immersed in magnetic fields above about 
0.1/zG and if the sources reside in magnetized clusters 
and super-clusters: For particles above the GZK cut- 
off this is because sources outside the Local Supercluster 
are beyond the GZK distance. On the other hand, sub- 
GZK particles are mainly confined in their local mag- 
netized environment and thus exhibit a much higher lo- 
cal over-density than their sources. Further, the sup- 
pressed flux of low energy particles leaving their envi- 
ronment is largely kept away from the observer if he is 
surrounded by significant magnetic fields 's^. Both ef- 
fects can be understood qualitatively by matching the 
flux j{E) in the unmagnetized region with the diffu- 
sive flux —D{E)\/n{E,r) in terms of the diffusion co- 
efficient D{E) and the density n{E, r) of particles of en- 
ergy E which shows that the density gradient always 
points to the source. More quantitatively, the shape of 
the large-distance component is demonstrated in Fig. 
which shows the observable flux resulting from an E~'^-* 
spectrum injected isotropically at a sphere with a ra- 
dius of 40 Mpc around the observer. Note that despite 
the smaller energy losses the sub-GZK particles arriving 
from outside the Local Supercluster are likely to have a 
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K„b,{ 40. < E:/B]eV < 9000.)= 57, £= 0.146 
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FIG. 6: Similar to Fig. 13 but for an observer in a much lower 
fleld region, ~ 8.2 x 10~^^ G. This corresponds to scenario 
2 in Tab. The overall likelihood significance is ~ 0.15 for 
n = 4 and < 10° in Eq. ©. 
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100 
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FIG. 7: The unnormalized energy spectrum observed at Earth 
resulting from a * isotropic proton flux injected at a 
sphere of 40 Mpc radius around the observer for scenario 1 
from Tab.0 averaged over solid angle and 10* computed tra- 
jectories. In contrast to the spectrum of the local component 
shown in Fig. 2] there is a clear tendency that this cosmolog- 
ical component can flt the flux neither at the highest nor at 
the lowest energies. 



spectrum even more strongly suppressed than in Fig. [3 
at low energies due to their containment in the source 
region. A significant contribution from sources at cosmo- 
logical distances to sub-GZK energies thus requires that 
neither these sources nor the observer are immersed in 
too strong magnetic fields and/or an injection spectrum 
considerably steeper than E~^ '^ to compensate for the 
systematic suppression of flux of lower energy particles. 

The confidence levels that can be obtained with this 
method for specific models of our local magnetic and 
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UHECR source neighborhood will greatly increase with 
the increase of data from future experiments. Full sky 
coverage alone will play an important role in this context 
as many scenarios predict large dipoles for the UHECR 
distribution. This is the case for basically all scenar- 
ios considered here, as demonstrated in Fig. |H1 Whereas 
current northern hemisphere data are consistent with sce- 
narios with Ns > 10 nearby sources at the ~ l.Ssigma 
level if the observer is surrounded by relatively strong 
fields i? ~ 0.1 /iG, a comparable or larger exposure in the 
southern hemisphere would be sufficient in these cases to 
find a dipole at several sigma confidence level, as demon- 
strated in Fig. |S| 



150 events with 40. < E/EeV < 9000. 



0.0100 



o 0.0010 




0.0001 



FIG. 8: Same as Fig. H but for comparison of the model 
predictions with an isotropic distribution [horizontal line, 
Ci ~ (47rA'')~^, see Eq. for the full-sky detector a la Auger 
discussed in the text, for N = 150 events observed above 40 
EeV. 



of Ci for the isotropic and the anisotropic distribution, 
respectively. Therefore, at ~ 10^^ eV an isotropic cosmo- 
logical flux about a factor 3 higher than the anisotropic 
flux originating within ~ 50 Mpc would be needed to ex- 
plain the isotropy observed by AGASA. For charged pri- 
maries this implies steep injection spectra and/or weak 
magnetic fields around observer and sources, as explained 
above. Without going into a more detailed analysis we 
remark that this will also require to decrease the flux 
contribution from nearby sources shown in Fig. 0] at the 
low energy end. As a consequence, the best fit injection 
spectrum for the local component will be slightly harder 
than the power law indices shown in Tab. ^ This is 
consistent with what is expected from shock acceleration 
theory |4|. 



1500 events with 10. < E/EeV < 9000. 



0.0010 



0.0001 




FIG. 9: Same as Fig. El but for iV 
above 10 EeV. 



1500 events observed 



Finally, the distributions of events down to 10"'^ eV also 
contain important information. Fig. |51 shows the multi- 
poles predicted by our standard scenario 1 in Tab.Qlthat 
full-sky experiment would observe for 1500 events de- 
tected above 10^^ eV. This corresponds to twice the num- 
ber of currently observed AGASA events and thus ap- 
proximately reflects the current exposure. A correspond- 
ing figure for the AGASA detector alone would look simi- 
lar. It is obvious that there is significant anisotropy even 
at I ~ 10, inconsistent with current AGASA observa- 
tions. On the other hand, cosmic variance becomes more 
important at these lower energies, and a possible signif- 
icant contribution from large-distance sources cannot be 
excluded if their magnetization is not too high, as dis- 
cussed above. It is easy to see from Eq. (0} that if a frac- 
tion fa of N events observed stems from an anisotropic, 
local contribution, whereas the fraction 1 — /a is cosmo- 
logical and completely isotropic, then 

ci - a, ((1 - ur + ^.n) , (7) 

where Ci^i = {AttN)~^ and Ci^a are the expectation values 



IV. CONCLUSIONS 

In the present work we performed UHECR propaga- 
tion simulations based on the distributions of magnetic 
field and baryon density obtained from a simulation of 
large scale structure formation. The magnetic field was 
simulated as a passive quantity and normalized at sim- 
ulation end in agreement with published measurements 
of Faraday rotation measures for groups and clusters of 
galaxies [23|. We considered finite numbers of discrete 
UHECR sources with equal total power and injection 
spectrum. Their positions were randomly selected with 
probability proportional to the baryon density. The ob- 
server was chosen within small groups of galaxies char- 
acterized by gas temperatures around a fraction of a 
keV, typical for our local environment. One chosen ob- 
server was found in a relatively high field region with 
B ~ 0.1/xG. For comparison, we also chose an observer 
situated in a small void, where the surrounding field is 
B ~ 10"" G. We found that good fits to the AGASA 
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data above 4 x 10 eV in the Northern hemisphere are 
only obtained for Ng > 10 sources and for observers sur- 
rounded by ~ O.l^G fields. Otherwise the predicted ar- 
rival direction distribution is either too anisotropic or 
produces too large auto-correlations at angles larger than 
a few degrees. The best fit case occurs for Ns ~ 100, sig- 
nificantly higher than in previous work 39] due to the 
more localized and more strongly structured magnetic 
fields considered here. 

For the required local source number density and con- 
tinuous power per source we find fi-sourco ^ 10~^ — 
10~^ /ig7Mpc~^, and Qsourco ^ 5 x IC'^ergs"^ respec- 
tively, the latter within about one order of magnitude 
uncertainty to both sides. This corresponds to an av- 
erage UHECR emissivity of quHECR = "-source Qsourcc ~ 

XO^^erg Mpc^'^ s~^ also with an uncertainty of roughly 
one order of magnitude, not larger, since it is fixed by 
the observed UHECR flux. 

Possible sources marginally consistent with these en- 
ergy requirements are radio galaxies. Their present en- 
ergy release of gig 5 x 10'^^ ergs~^ Mpc"'^ (/powcr/10) 
|49j is roughly what is required in order to produce a suf- 
ficient flux of UHECR, assuming that the injection power 
law is flat (oc [^|^. The parameter /power de- 

scribes the ratio of the total power of the radio galaxy to 
the equipartition estimate based on its radio luminosity, 
and it enters the used radio luminosity-jet-power relation 
of We expect /power 10 within an order of mag- 
nitude. In the estimate of the radiogalaxy power the 
observed radio luminosity function j 52 | was integrated 
only for sources with a 2.7 GHz luminosity of more than 
Lmin = 2 X 10^^ Watt Hz~^ , since they correspond to a lu- 
minosity of lO^'^ergs"^ (/power/10). This would provide 
the required UHECR luminosity per source of Qsourcc ^ 
5 X 10'^^ergs~V/nowBr /tihf.cr/O.S), using the optimistic 
assumption of |^ that /uhecr = 3 — 10% of the 
radio galaxy power is converted into UHECRs. The im- 
plied number density of these radio galaxies is nsourcc ^ 
10^"* Mpc^^ and, therefore, is only marginally consistent 
with the required n-souree ~ 10"'' — 10~^ /i|yMpc~'^. Since 
the number density increases strongly with decreasing 
imin this requirement can possibly be fulfilled by allow- 
ing for a larger number of less powerful UHECR sources. 
This implies that basically every radio galaxy has to be 
an efficient UHECR source, not only the most powerful 
ones. Since many of the weaker radio galaxies do not 
exhibit a hot spot, which is assumed to be the UHECR 
acceleration site in the scenario of 50, 51], their efficiency 
in producing UHECR might be largely reduced. This is 
a potential serious problem for this scenario, since low- 
ering /uhecr by several orders of magnitude can not 
be fully compensated by assuming a higher radio galaxy 
jet-power, because /power > 100 does not seem to be con- 



sistent with observations of radio galaxies j53] . 

To conclude, radio galaxies can be the sources of UHE- 
CRs if even weak radio galaxies are efficient particle ac- 
celerators to the highest energies, otherwise they have 
serious problems to reproduce the smooth UHECR ar- 
rival direction distribution. 

We also found that consistency with the isotropy ob- 
served by AGASA down to 10^^ eV requires the exis- 
tence of an isotropic component with a flux about a fac- 
tor 3 larger than the local component. This isotropic 
component would presumably be of cosmological ori- 
gin and thus would not contribute significantly above 
4 X 10^^ eV due to the GZK effect, consistent with the 
fact that at these energies we find local scenarios consis- 
tent with all data. The resulting best fit injection spec- 
trum for the local component is £'"(2.2-2.4)^ contrast, 
for the charged primaries of the cosmological compo- 
nent to dominate around lO'^^ eV steep injection spectra 
and/or weak magnetic fields around observer and sources 
would be required. These two conflicting requirements 
provide a strong argument against the hypothesis of lo- 
cal astrophysical sources of UHECRs above ~ 10^^ eV 
in a strongly magnetized and structured intergalactic 
medium. 

Finally, we have also demonstrated that already a 
modest increase in data together with full-sky cover- 
age will allow to put considerably stronger constraints 
on UHECR source and magnetic field scenarios than 
presently possible. In particular, our local scenarios pre- 
dict the emergence of significant dipoles and quadrupoles 
above 4 x 10^^ eV. 

Modeling our cosmic neighborhood and simulating 
UHECR propagation in this environment will therefore 
become more and more important in the coming years. 
This will also have to include the effects of the Galactic 
magnetic field and an extension to a possible heavy com- 
ponent of nuclei. For first steps in this direction see, e.g. 
Refs. and Ref. [H, respectively. 
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